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Abstract 

This letter announces and summarizes results obtained in [5] and considers several 
natural extensions. The aforementioned paper proposes a procedure to reconstruct coef- 
ficients in a second-order, scalar, elliptic equation from knowledge of a sufficiently large 
number of its solutions. We present this derivation and extend it to show which parameters 

VO ' may or may not be reconstructed for several hybrid (also called coupled physics) imaging 

modalities including photo-acoustic tomography, thcrmo-acoustic tomography, transient 

H , ' elastography, and magnetic resonance elastography. Stability estimates are also proposed. 

<' 

J3 \ 1 Introduction 
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a _ 

Consider a general second-order, linear elliptic equation of the form 

V • aVuj + b ■ Vuj + cuj =0, x £ X, Uj = fj x G dX, 1 < j < J, (1) 

for X a smooth open domain in M n , with n spatial dimension, and (a, b, c) possibly complex- 
valued, symmetric second-order tensor, vector field, and scalar coefficient, respectively. We 
assume that a is elliptic, the real part of a is coercive and bounded, and c is such that the above 
equation admits a unique solution. We also assume that (a, b, c, V • a) are of class C°' a for 
^r) \ some q > 0. We construct J solutions of the above equation for different boundary conditions. 

Several recent hybrid inverse problems aim to reconstruct the unknown coefficients (a, b, c) 
from knowledge of internal functionals of the coefficients and of the elliptic solutions (uj)i<j<j- 
Concretely, we assume knowledge of the following functionals 

k>( ■ Hj(x) = d(x)uj(x), x G X, (2) 

j^ I with d{x) a scalar coefficient that is a priori also unknown. 

What may be reconstructed from (a, b, c) from knowledge of {Hj)\<j<j in the setting d = 1 
is analyzed in [8]. We present the reconstruction procedure of the aforementioned reference 

IT. ,, . 

in section [2j Such a reconstruction is based on the availability of ratios of solutions j^- = ^- . 
This preliminary step is then used in section [3] to show that (a, b, c, d) can be reconstructed up 
to explicit obstructions that take the form of gauge transformations. We also provide stability 
estimates for the reconstructions. 

The inverse problems with internal functionals of the form ([2]) are part of a larger class 
referred to as hybrid inverse problems or coupled-physics inverse problems. For recent results 
and reviews on this topic, we refer the reader to pQ HI QUI [321 CEH CC2 [IH HU [El [TT] and their 
references. 
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2 Reconstruction procedure 

For the rest of the paper, we assume the existence of u\ 7^ on X . We refer to [8] for conditions 
on /1 that ensure such a property, either globally, in favorable cases, or at least locally. This 
allows us to define the known quantities 

^- = %i = ^±i, l<j<J-l. (3) 

Using the notation A : B = Tr(AB) for symmetric matrices A and B, we find that 



h 



a T : V®V- + T -Vv j = O x e X, Vj = ^±± x € dX, (4) 



where for an arbitrary complex-valued non-vanishing function r(x) on X, we have 

a T = ru-^a, fi T = Tu-^b + tV ■ aiii. (5) 

Note that the equation Q is invariant by multiplication by a non-vanishing scalar coefficient 
so that (a r ,/3 r ) may at best be reconstructed up to a multiplicative scalar coefficient. The 
result in [8] shows that this is the only obstruction to the reconstruction of (a T ,0 T ). 

More precisely, let us assume that (V?;i, . . . , Vu n ) form a basis of M n for all x £ X. We 
distinguish the case a scalar from the case a a second-order tensor. When a is scalar and 
J = n + 1, then (a T ,0 T ) are reconstructed up to the multiplicative scalar r. This is equivalent 
to saying that a~ l b is uniquely reconstructed. Indeed, we have 

b 

Av j + — • Vvj = Avj H Vvj = 

J a T a 

so that, defining Hj,- = Vfj • Vvj and H l i the coefficients of iJ _ , we have 

aT x b = H ij {a~ l b ■ V Vj )Vvi = -H ij A Vj Vvi, (6) 

where we have used the convention of summation over repeated indices and the fact that for 
any vector F, we have F = H 13 F ■ VvjVv j. 
When a is tensor-valued, we need 

J = I n := -n(n + 3) = n + 1 + M n , M n = -n(n + 1) - 1. (7) 

For 1 < j < I n — 1 and 1 < m < M n , let us define the coefficients QJ 1 such that 

In-l In-1 

^efVvj = and the symmetric matrices M m = ^ O^V^Vj, (8) 

such that (M m )i< m <M n form a free family of symmetric matrices. Sufficient conditions are 
presented in [8] to guaranty that (Vu,)i<j< n and (M m )i< m <M n are free families for the choice 
gm _ -Hi k Vv m +n ■ Vufc for 1 < j < n, 9™ = 1 for j = n + m and 0™ = otherwise. The 
above construction allows us to obtain the following constraints: 

a T : M m = 0, 1 < m < M n . (9) 

This implies that a T = M°, where (M )* is a matrix in the one-dimensional orthogonal 
complement to (M m )i< m <M n for the inner product for symmetric matrices (^4, B) = Ti(AB*). 
Thus a T is reconstructed up to a multiplicative scalar coefficient. From (j4|), we deduce that 

T = -H ij a T : V® 2 VjVvi. (10) 



Note that the above is nothing but (|6j) when a is a scalar coefficient. 

This shows that (a T ,f3 T ) are uniquely reconstructed up to the multiplicative coefficient r. 
Note that additional information of the form H^ = du^ for u^ solution of ([1]) with u/, = fk on 
dX does not provide any new information. Indeed, -^ is a solution of the elliptic equation 
(j4]) with known boundary condition ^ on dX. 

3 Reconstruction of (a, b, c, d) up to gauge transforms 

Reconstruction up to gauge transforms. The above derivation shows that all that can be 
extracted from an arbitrary large number of functionals of the form H^ = duk is (a T , (3 T , Hi) 
augmented with the equation for u\. Let us decompose a = B 2 d for a a matrix with determi- 
nant equal to 1. We assume here to simplify that such a decomposition is valid globally on X 
(which is obvious in the case where a is real- valued and positive-definite). Since a T = tu\B 2 cl 
is known, we deduce that a is known. We compute 

aa~ 1 (/3 r - V • a T ) = — ^ — (Vlnr) • a. 

Moreover, defining v = Bu\ = — V- , we find that 

Av _ V • dVB c 
~ ~ B + B*' 

Thus, we obtain after elimination of r and u\ that knowledge of (a T , f3 T ,H±) and the equation 
for u\ is equivalent to knowledge of 

b , B A^- V-aVB c \ 

>^ + 2aVln 7 ,^ = ^^ + ^). (11) 

d 

No additional information may be extracted from functionals of the form Hk = duk since 
knowledge of the above coefficients uniquely determines the functionals H^- 

The dimension of the unknown coefficients in (jlip is — j — — l + n + l = I n = ^n(n + 3), 
which is the number of functionals used to reconstruct them. The dimension of (a, b, c, d) is 
—^2 — - +71+1+1 = in+2. There are therefore two gauge parameters that remain undetermined. 
Moreover, (a, b, c, d) are reconstructed up to any transformation that leave the coefficients in 
dill) invariant. 



Applications to medical imaging modalities. In the setting of Transient Elastography 
and Magnetic Resonance Elastography, we may assume that d is known (and equal to 1) and 
that 6 = 0. We thus obtain a (redundant) transport equation for B (or equivalently for the 
gauge r) and then an explicit expression for c. Therefore, (a, c) is uniquely reconstructed. 
More generally, when V • (a -1 b) is known, we obtain an elliptic equation for B or equivalently 
for r. Then (a,a~ 1 b,c) is uniquely reconstructed. 

In the setting of quantitative photo-acoustic tomography (QPAT), we may assume that 
b = and that d = Tc. We again obtain that -^ = j^ is known, and hence q = — ^, is known. 
The reconstruction of (B, c, T) is unique up to any transformation that leaves (-^, ^ — h -§z) 
invariant. When T is known, then (B,c) are uniquely reconstructed [HEJE]. 

A similar result may be obtained in the imaging modality called quantitative thermo- 
acoustic tomography (QTAT), where d = T(^sc)ul; see [21 [3j [7] for a derivation of such a 
model for Hj = T(^sc)uju*. Assuming again that b = 0, or more generally that V • (a" 1 ^) 
is known so that r, or equivalently -r is known, then v = Bu± is known. In this setting, we 



thus find that (B,c,T) are reconstructed up to any transform that leaves (^f, Va ^ B + -^-) 
invariant. Note that when a is real-valued, then V is uniquely reconstructed and (B,c) are 
reconstructed up to a transform that leaves Vl ^ 7 - B + -j^ invariant [3] . 

Note that, more generally, one condition on the field b is sufficient to uniquely reconstruct 
the gauge r or equivalently -j. Indeed, we observe that the second known quantity in (fTTj) is 
equivalent to knowledge of a~ l b + 2Vln -j. Thus, knowledge of one component of a~ 1 b, or of 
V • a~ 1 b, for instance, again provides an equation that allows us to uniquely reconstruct -y 
and, hence, a~ l b. In such a setting, q = —p with v = Bu\ = — ^— is known and (B,c,d) can 
then be reconstructed up to any transform that leaves (-^-, — ^ \- -^j) invariant. 

4 Sufficient conditions and stability estimates 

Sufficient conditions. The results of the preceding section exactly characterize which co- 
efficients in (a,b,c,d) can be reconstructed. Such reconstructions hinge on the solutions (uj) 
to be sufficiently independent. More precisely, we assume that u\ 7^ on X, (Vvj)i<j< n is a 
basis of M n at every point x G X, and that the matrices M m are linearly independent on X . 
In some situations, for instance when complex geometric optics (CGO) solutions can be 
constructed, the above conditions are shown to hold for an open set of well-chosen boundary 
condition (fj)i<j<i n [SJ- However, in the general situation where a is possibly complex-valued 
and anisotropic, such CGO solutions are not available. The linear independences mentioned 
above can be shown to hold locally on subdomains on X. More precisely, it is shown in 
[8] that for a finite covering U^^X^ of X, then for an open set of boundary conditions 
(fj)i<j<j=Kxi n , we can construct a non-vanishing solution Uk,i on Xk, linearly independent 
gradients C^^ hL )2<j<n+i and linearly independent matrices M k,m as constructed in ©. 

Stability estimates. The procedure leading to the reconstruction of (|lip is explicit and 
allows one to estimate how errors in the functionals (Hj) propagate into errors in the re- 
constructed coefficients. Let us assume that d is known and smooth and that 6 = for 
concreteness. Similar results can be obtained in more general cases. We observe that the 
construction of the matrices M m involve taking two derivatives of the functionals Hj. The 
reconstruction of o therefore involves differentiating (Hj) twice. 

When b = 0, we observe that the reconstruction of V-B or equivalently Vr from the (second) 
vector field in (jlip also involves differentiating (Hj) twice. Once (d,B) are known, then (llip 
provides a formula for c. However, some simplifications occur. From ([5]), we observe that V • a 
is reconstructed from differentiating (Hj) twice (and not thrice). Then with u\ known since d 
is known, we reconstruct c directly from (pQ) with again a loss of two derivatives. This yields 
the result 

\\(a,c, V • a) - (a,c,V ■ a)\\ c o, a + \\B - B\\ c i, a < C\\(Hj - H j ) 1 <j< J \\ C 2, a , 

for some positive constant C, where Hj is constructed as Hj in ([2]) with the coefficients (a, b, c) 
in (pQ) replaced by (a, 0, c). Similar stability estimates may be obtained in the more general 
case with d and b unknown; see [8] for additional results. 
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